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ABSTRACT 

In this paper the problem of an elastic half plane containing a crack 
and stiffened by a cover plate is considered. First, the asymptotic nature 
of the stress state in the half plane around an end point of the stiffener 
is studied in order to determine the likely orientation of a possible frac- 
ture initiation and growth. The problem is then formulated for an arbitrary 
oriented radial crack in terms of a system of singular integral equations. 
For an internal crack and for an edge crack, the problem is solved and 
the stress intensity factors at the crack tips and the interface stress are 
calculated. The case of a cracked half plane for two symmetrically loca- 
ted cover plates is then considered. From a fracture viewpoint, the case 
of two stiffeners appears to be more severe than that of a single stiffener. 


1. Introduction 

In the relatively recent past, there has been considerable interest 
in the analysis of "cover plates" as a problem in structural mechanics. 
This is primarily due to the fact that the structural components with a 
variety of bonded or welded stiffeners and the solid state devices con- 
taining elastic wafers fuse-bonded to elastic substrates may be approxi- 
mated by a cover plate bonded to an elastic solid. In most cases, since 
the stiffener is relatively very thin compared to the remaining dimensions 
of the structure, it is approximated by an elastic "membrane" neglecting 
the normal stress along the interface. Typical examples for such studies 
may be found in [1-3]. The primary interest in these and similar studies 
has been in the evaluation of the contact shear along the interface. The 
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results given in [3] for various elastic and inextensible cover plate 

geometries show that, at the end points of the stiffener, not only the con- 

- 1/2 

tact shear stress but also the stress state in the substrate has a r ' 
singularity. This suggests that such points of stress singularity are 
locations of potential failure initiation. Furthermore, if the bond is 
sufficiently strong, the most likely failure mechanism would be the initia- 
tion and propagation of a crack in the substrate along the weakest cleavage 
plane emanating from the singular point. To study the related failure 
problem, one needs to determine the weakest cleavage plane in the elastic 
half plane and to solve the corresponding crack problem by placing a crack 
of arbitrary length along this plane. 

A problem similar to that described above was considered in a recent 
paper [4], where it was assumed that the stiffener is perfectly rigid and 
the crack is perpendicular to the boundary. Aside from the assumption 
regarding the stiffness of the cover plate which may not be very realistic 
for most practical applications, in [4] it is found that the nature of the 
stress singularities at both ends of the cover plate are identical to 
that of stiffened uncracked elastic half plane where the power of singu- 
larity is complex. This ‘is clearly incorrect as the zero traction-zero 
displacement mixed boundary conditions which prevails in the half plane 
near the end points are for a wedge of angle ir at one end and ir/2 at 
the other At the crack end of the stiffener, that is, for the 90- 
degree wedge, the power of the singularity is less than one half and is 
real [5]. 

The general formulation given in this paper may be used to reduce a 
number of crack-stiffener or crack-contact problems (including that dis- 
cussed in [4]) to a system of integral equations. Since the kernels in 
these integral equations are known in closed form and are relatively 
simple functions, quite accurate solutions to the problem may be obtained 
with a minimal computational effort. 

This is due to "rounding" the corners by approximating the mapping 
function by a rational function used in [4] for mapping the cracked 
half plane into a circle. 
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2. Formulation of the Problem 

The problem under consideration is that of an elastic half plane 
(-00 < X < ~, y < 0) which contains an arbitrarily oriented crack and is 
subjected to a given set of external loads. Along the boundary y = 0 
the plane may contain a bonded stiffener or may be loaded through a stamp. 
Perhaps the simplest way to solve the problem would be to reduce it to a 
system of integral equations in which the crack surface displacements and 
the contact stresses on the boundary are the unknown functions. To derive 
the integral equations, in addition to the solution of the problem for the 
half plane under the given applied loads but without the crack and the 
stiffener, it is sufficient to obtain the solution for a dislocation in 
the plane and a concentrated force on the’ boundary. 

Consider a half plane (-<» < x < <», -<» < y < 0) with the elastic con- 
stants < 2 * ^2 (<2 ~ ^ ~ plane strain and <2 = (3-v2)/(l+v2) 

plane stress) which contains a dislocation at the point (x-j, y-j) having 
a Burger's vector S. Let the components of the Burger's vector be b = 

X 

f-j and by = f 2 « Referring to [6] the stress state in the half plane due 
to f-j and f 2 may be expressed as 

2^2 

oixx(x,y) = [K^-|(x,y,x.|,y.|)f-] + K-j2(x,y,Xpy^)f2], 

2p2 

oixy(x>y) = Tr(i+i<2) t*^21^^’^’^l’^l^^l ^ 

2^2 

0lyy(x,y) = [K3l(x,y,Xi,yi)fi + K32(x,y,x-| ,y^)f2l , 


(l.a-c) 


where the functions are given in Appendix A. 

The second basic solution needed to formulate the problem is that 
of a half plane under concentrated forces f^ and f^ acting on the boundary 
at (x = Xq, y=0) in x and y directions, respectively. This solution is 
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given by 


= - :jjpr C(x-Xq) + (x-XQ)^y f^], 

*^ 2 xy^^’'^^ irr‘* ^3 ^ 4 ]* 

»2yy(x.y) = - ^ [y^(x-x^)f3 + y^f^]. (2a.c) 

where 

= (x-x„)2 * / . (3) 

For this loading condition, the displacement derivatives on the boun- 
dary y=0 may be expressed as 


4uo a 1 fo 

Tfi^BT '7^ 

d 0 


K 9-1 

+ — = f 

<2+1 4 


5{x-Xq) 


'^'>z a 

I+ 1 C 2 8X 


V2(x,0) = - 


“ 1 

— = f 

<2+1 3 


6(X-X ) + - — ^ . 
O'' TT X^-X 


(4a, b) 


The third problem is the determination of the stress state in the 
half plane having no crack and no stiffener and subjected to the actual 
applied loads. This stress state will be designated by Oj^-{j{x,y), 

(i»j = x,y). For example, in a plane under uniform tension p^^ in x- 
di recti on, we have 


^^axx^^’"^^ ^ 0 ’ '^axy^^*'^^ '^ayy^^*'^^ ~ (5a-c) 

Consider now the problem described in Figure 1. If the cover plate 
is approximated by a membrane, 0 ^^ = 0 and the equilibrium condition in 
x-direction gives 
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X 


* '’l 


4 


8^1 du-j 

1+K^ ax 


h = 


-2a 


( 6 ) 


where y-j and k-j are the elastic constants, h the thickness, u-i the dis- 
placement, 2a the length of the stiffener, P-j and ?2 the additional forces 
which may be acting on the stiffener, and f 3 (xQ) = the shear stress 
acting on the contact region (Figure 1). Note that 
0 

I f3{Xo)'‘Xo ' '’2 ■ ’’l • (7) 

-2a 

If we now assume that the half plane contains a crack along the 
line L, the Burger's vector ^(f-| ,f 2 ) would be a continuously distributed 
function with the coordinates x-j and y^ on L. Clearly, f^(xpy-j), 

unknown functions of the problem which may 
be determined from the traction boundary conditions on the crack surface 
and the displacement continuity along the contact area. Anticipating 
the crack initiation at an end region of the stiffener, consider the 
specific radial crack geometry shown in Figure 1. The boundary condi- 
tions on the crack surface may be expressed as 

cTnn(s) = a^j^cos^e + a^ySin^e + o^ySin2e = 0, c < s < d, 

asn(s) = {a^^-Oyy)sinQ cose - Cj^yCos2e = 0, c < s < d, (8a, b) 

where (see (1), (2) and (5)) 

®ij ' °lij + “21j ^ '’alj- <’’j> ' 

On the line of the crack note that 

X = s sine , y = -s cose , 


(9) 


( 10 ) 
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x<| = t sine, y-f =-t cose, (c < t < d). 


( 11 ) 


Now, observing that f-|, f2 and are distributed functions, substituting 
from (1), (2), (5), and (9-10) into (8) we obtain the following two 
integral equations: 

d d 0 

I k^^(s,t)f-j(t)dt + I j 

c c -2a 

= - Pq cos^e, (c < s < d), 
d d 0 

I k2-j(s,t)f^(t)dt + j k22(s,t)f2(t)dt + j k23(s,XQ)f3(xQ)dx^ 
c c -2a 


= - Pq sine cose, (c < s < d) , (12a, b) 

where the kernels k.., (i=l,2; j=l,2,3) are given in Appendix B. Since 

* J 

the crack is an embedded crack, from the condition of single-valuedness 
of displacements it follows that 



d 

f2(t)dt = 0 . 
c 


(13a, b) 


The third integral equation is obtained by expressing the condition 
of continuity of 3u/3x along the interface y=0, -2a < x < 0, namely 


The strain 3U^/3x is given by (6). 3U2/3x is obtained by adding the 

strains given by the stress states (1), (2) (which is given by (4a)), 



and (5), and again keeping in mind that fg, and are distributed 
functions. Thus, for the crack geometry and the applied loads given in 
Figure 1 we obtain 



d 0 

k32(x,t)f2(t)dt + I k33^^*^o^'‘'3^^o^^^o 
c -2a 


l+<2 

Pq SinTj' *"1 


(-2a < X < 0), 


(15) 


where the kernels k 3 j, (j=l,2,3) are also given in Appendix B. 


3. The Uncracked Plane 

To determine the direction of crack initiation in the plane, the 
asymptotic stress state near the end points of the stiffener needs to be 
analyzed. In the uncracked plane, the external loads are the actual 
applied loads (e.g., a„„(«,y) = p.) and the contact stresses on the 
stiffener-half plane interface. If the stiffener is approximated by a 
membrane, the contact stress ayy(x,0) = f^Cx) is zero and referring to 
the insert in Figure 2, for P-|=P2=0, o^y(x,0) = f 3 (x) is determined 
from (see (15)) 



x„-x 


1 

^0 2h y-| (1+K2) 


X 

0 


subject to 
b 

0 


Pq/2, (0 < X < b), 
(16) 

( 17 ) 


Defining 
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( 18 ) 


bp2Cl+K-|) 

• X ' (l+t)b/2, x^ = (l+x)b/2, fjCXg) = f(r). X = • 

ec|U3t1ons (16) and (17) may be expressed as 


1 

1 ■ 

IT 

-1 


^(0- 


dx 



- Pq/ 2, (-1 < t < 1) , 


1 

» 

f(T)dx = 0 

-1 


(19a, b) 


which are solved numerically. 

Once f3(xQ) is determined, by substituting from (2) into the trans- 
formation fomulas (8) and integrating in x^, one may easily obtain 
the cleavage stress along n=0 as follows: 


J_ 

Po 



1 

i ■ 

IT 

-1 


{[r sine + (l-r)/2]^ 


cos^e + [r sine 


+ (l-T)/2]r2 sin^e cos^e - 2r[r sine 
+ (l-x)/2]^sine cos^eHr^cos^e + [r sine 
+ (l-x)/2]^}“^f(x)dx + cos^e , (0 < r < «), 


( 20 ) 


where r = s/b. Also, on the boundary y=0 the stress component may 
be obtained from (2a) as follows: 


— a (x 0) = 1 + - f dx , X = b(t+l)/2, (-» < t < “). (21) 

P XX' ’ ' IT J T-t 
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Note that the solution of (19a) is of the form 


1/2 

f(x) = F(t)/(1-t2) (-1 < T < 1). (22) 

F(t) being a bounded function. Thus, after determining F(t) from (19), 
an„(r,e) and ay„(x,0) may be obtained from (20) and (21) by using the 
standard Gaussian integration formulas [7]. Figure 2 shows the variation 
of the cleavage stress with the angle 9 for a fixed value of r=s/b. 

The value of r seems to affect primarily the magnitude rather than the 
angular variation of The angle at which is maximum is approxi- 

mately 3 degrees. By changing x it was also observed that this angle 
does not seem to vary significantly with material constants. Therefore, 
for the remainder of this study, it will be assumed that the weak cleavage 
plane in the elastic half space is e = 0. This way the integral equations 
and the subsequent asymptotic analysis are simplified quite considerably 
without significantly altering the qualitative behavior of the results. 
Along the plane e = 0, the variation of the cleavage stress aj^j^(s,0) = 
Ovv(b,y) is shown in Figure 3. From (20) it can be shown that a has 
a singularity of the form s*'^ in the neighborhood of the end points of 
the stiffener. Thus, becomes unbounded as s 0. 

By substituting from (22) into (21) and by using the properties of 
the Cauchy-type integrals or that of the Chebishev polynomials, it may 
easily be shown that [7] the function defined by (21) is 

bounded in -1 < t < 1, (or in 0 < x < b) and has a square-root singularity 
at t = + 1 (or at x=0, x=b) for |t| > 1. That is, as may be seen from 
Figure 4, cTy„(x,0) is discontinuous and is unbounded at the end points 
of the stiffener. 

4. The Crack-Stiffener Problem-The Internal Crack 

For uniform tension p^, P-| = P 2 “ 0 = 0 from (12) and (15) 

the integral equations of the problem shown in Figure 1 may be obtained 
as follows: 
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] fi(t)dt 


iiTT+i^ 



0 



J 1_ _ 2t , 

t-s ~ t+s (t+s)^ 


x.3f (x )dx. 

0 3' 0' 0 _ 

(Xo^+S^)’^ “ ■ 



4t2 

Tt+?F 


(c < s 


< d) , 


( 23 ) 


^^2 f r 1 4. s2-t^+4tsi 

ir(l+K«) J ^ t-S (t+S -* 

C 


f2(t)dt - I 



(Xq2+s^P 


9 


(c < s < d) 


(24) 





, '-^2 t° 

%,„J_2a x„-x ‘”'0 


1+K 

Shir 





1+IC2 

8i^'’o • 


(-2a < X < 0) , 


(25) 


The solution of (23-25) is to be obtained under the following condi- 
tions: 



d 0 

f2(t)dt = 0» j ■•'3(^o^'^^o "" ° * 

c -2a 


(26a-c) 


For oO the integral equations (23-25) have simple Cauchy-type 
kernels. Consequently, the functions f. , (i=l,2,3) have square root 
singularities and are of the form shown by (22). These equations may 
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easily be solved by using a Gaussian quadrature formula [7]. From the 

viewpoint of fracture analysis, the quantities which are of primary 

interest are the stress intensity factors which may be obtained as 
(*) 

follows' . 


k^(c) = lW2(c-s)aj^j^(s,0), k2(c) = lim*^(c-s)axy(s,0) , (27a, b) 


s-k: 


s->c 


k-|(d) = linn/2(s-d)a^^(s,0), k2(d) = lim»^2ts-d)a^y(s,0) 


s-kI 


S-Ki 


k2(-2a) = lim/2(2a+x)a (x,0), kp(0) = lim/^a^^(x,0). (29a, b) 


x-^-2a 


x^- 


(s,0) , 

(27a,b) 

(s,0). 

(28a, b) 

(x,0). 

(29a, b) 


5. The Edge Crack 

The integral equations (23-25) are valid for all cracks perpendicular 
to the boundary, including the physically important case of an edge 
crack for which c=0. It may be observed that for c=0 the kernels of the 
integral equations are of the generalized Cauchy type; that is, in addi- 
tion to simple Cauchy singularities, the kernels contain terms which 
become unbounded at the end point x=0=s. Consequently, at the end point 
x=0=s the functions f^. , (i=l,2,3) v/ould not have the standard square- 
root singularity. 

First, we note that for a sectionally holomorphic function F(z) 
defined by 


■pn 

' 'Note that for the geometry under consideration 

^l^y^ [u2(+0,y)-U2(-0,y)] and fg(y) = [v2(+0,y)-V2(-0,y)], 

s=-y, <^ns~"Sy’ transforming the coordinates the notation f^(s) = 

f,*Cy(s)], (i=l,2) is used. 
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( 30 ) 


b 



f (t)dt 
t-z 


where 


f(t) 


(t-a)“(b-t)» 


(-1 < Re(a,e) < 1), 


(31) 


one may express the following asymptotic relation [7,8]: 


F(z) = 


»(a) e^°^ 1 

(b-a)^ (z-a)“ 


j»(b) 1 1 

(b_a)« sin$Tr 


+ Fq(z) , (32) 


where ((i(t) is a bounded function and is either bounded or has singular- 
ities of order lower than that of F(z) [8]. Also note that 


’ J (t-z)2 

a 


dt = ^ FU) 



f(t) 

(t-z)^ 


dt 


dz2 


F(z) , 


(33) 


By observing that for a < x < b F(z) is holonorphic at z=2a-x, one 
could write 


1 

TT 


I 


f(t)dt 1 ^ 
(t+xr2a)'^‘^^ ^ dz" 


F(2a-x) 


(a < X < b) . 


(34) 


Thus, for c=0 the asymptotic values of all the integrals in (23-25) may 
be expressed in terms of the corresponding holomorphic functions by using 
the general relations (32) and (33) and specific expressions of the form 
(34). For this, it is sufficient to expand the kernels into simple 
fractions. For example. 
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Xt2 _ i _ 1 1 1 ^ X 1 X 1 

(t2+^2)2 4 tMx 4 t-lx 4 4 (J .,^,2 - 


( 35 ) 


for which the related function F 2 (z) would be holomorphic at z ~ + ix if 
-2a < X < 0 (see (25)). 

We now let the unknown functions f^., (i=l,2,3) in (23-25) be 
defined by 


2^2 , ^ <f>i(t) 


(36) 


2Un <{)«(t) 


<l>3(t) 

^3(^) " / o \a ’ < t < 0), 0 < Re(a,3,Y,u) < 1 . (38) 

l-tj (t+2a) 


Also let 


Fj(z) 


2^2 1 f 

~ l+Kj, 7 J “t-z (j=l»2) , 


f<(t) 


(39) 


F,(z) 


1 fO ^3^*) 


- 1 f 

" IT J t 
-2a 


dt . 


(40) 


By substituting from (36-38) into (23-25) and by using the function- 
theoretic method described above, the following characteristic equations 
may be obtained to determine the constants a, b, y> and oi: 

‘fr-jCd) <l>o(d) <f»o(-2a) 

— s — cotTTO) = 0, — cotiTY = 0, — cotiTct = 0, (41a-c) 

d^ d^ (2a)^ 
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1 * ^i(o) o^^oCo) 

BnJs [0-«S’'6-‘l6+262) — + (B-2)cos f ^ ] = 0 , (42) 

*1 ^oCO) Q 4^o(0) 

^^[(cos.6-H4s-2b 2) -^+ (B-l)slnf--^] = 0. (43) 

d' (2a) 

0 . 

(44) 

At the end points, the functions <j>^.(t), (i=l,2,3) are bounded and 
nonzero. Hence, equations (41) give the following expected results cor- 
responding to square-root singularities: 

to = 1/2 , Y = 1/2 , a = 1/2 . (45) 

Since it is assumed that <|)^.(0) f 0, (i=l,2,3), the coefficient determinant 
of the linear homogeneous system (42-44) must vanish, giving the fourth 
characteristic equation to determine B as follows: 

[costtB - 1 + 4b - 2b^]^ = 0. C46) 

sin^TrB 

It may be observed that (46) has no root for which 0 < Re(B) < 1. 

For the sectional ly holomorphic function F(z) given by (30), if we now 
assume that the density function f(t) is bounded at t=b and let 

f(t) = , (47) 

(t-a)"* 

the asymptotic expression of F(z) near the point z=b becomes [8] 


_J 

SinirB 


[-2BCOS ^ ^ + 2(1-B)sin ^ — cosirB — ] = 


<i>i(0) 


ttB '"2 


4*9 ( 0 ) 


<^o(0) 


d^ 


(2a)“ 
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F(z) = ^ log (z-b) + F^(z) . 


(48) 


where, near and at z=b, is bounded. With the behavior at the end 
points known, the density functions f^. may now be expressed as 


2^2 , , 


2^2 

<2+1 


f2(t) = 


(d-t)** 


» f 3 (t) = 


(t+2a)^ 


( 49 a-c) 


The functions G-j and G2 at t=d and G2 at t=- 2 a are bounded and nonzero. 
However, their behavior at t =0 is as yet unknown. By using ( 48 ) and 
substituting from ( 49 ) into the integral equations ( 23 ) and ( 25 ), we 
obtain 

-L 63(0) log X < - , - 63(0) log s < » . ( 50 a,b) 

irVZa TTI^^ 

In the third integral equation ( 24 ) the coefficients of the logarithmic 
terms cancel out. At the end point, we have x = 0 - and s = O’*". Hence, 
from ( 50 ) it follows that 


63(0) = 0 . ( 51 ) 

Therefore, the characteristic equation ( 46 ) found by assuming that 
<1)3(0) 0 (which means that f3(0) 0 or G3(0) ^ 0) would not be valid. 

Going back now to the system of equations ( 42 - 44 ), if we let 413(0) = 0, 
and assume that 4 .^( 0 ) 0 and 4>2(0) ^ 0 ( 42 ) and ( 43 ) gives 

(1 - cosxe - 4e + 2$2) = 0 , ( 52 ) 

and from ( 44 ) we obtain 
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4'iCO) ~ 2 ’ 


(53) 


It should be pointed out that (52) is the characteristic equation 
corresponding to a 90-degree elastic wedge for which the tractions in 
the neighborhood of the apex are zero. In the problem under considera- 
tion, since a (x,0) = 0 and (^JQ) = 0 implies that a^,,C-0>0) = 0, this 
result is expected' . The relation (53) indicates that even though 
bounded and nonzero, at (x=0, y=0) the displacement derivatives 
3(u2'*’-u£)/3y and 3(v2‘’’-V2")/3y are not independent. In passing one may 
also remark that at y=0, if one forces the crack to close (i.e., if 
(j)-j(O) = 0 and (}>2{0) = 0) and lets ({> 2 ( 0 ) f 0, from (44) one may easily 
recover the standard characteristic equation, namely cotirB = 0. 


6. Symmetric Cover Plates. 

• The formulation of the stiffener-crack problem described in the pre- 
vious sections may be applied to a cracked half plane containing any 
number of stiffeners without any difficulty. In particular, the problem 
is considerably simplified if there are two identical cover plates 
located symmetrically with respect to the x=0 plane and if the crack is 
oriented along the plane of symmetry. In this case, fg(t) = 0 and the 
system of integral equations (23-25) reduces to 


5 d 

^^2 f r_ J 1_ 

Tr(l+tC2) J ^ t-S t+S 


(t+s)2 (t+S)3 ’ 


__ 

' 'It should be anphasized that the absence of stress singularity at 
the apex of a 90-degree elastic wedge stiffened on one side is 
due to the membrane assumption made for the stiffener. If the 
stiffener has a finite thickness, then the apex is a point of stress 
singularity and the asymptotic behavior of the stress state around 
this point is similar to that of two bonded dissimilar 90-degree 
elastic wedges (see, for example, [5]). 
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+ 


-b 

i f t3 

' -2a tt2+s2)^ 


fjCtjdt = -p^ 


(c < s < d) , 


C54) 


I i L 'A • 

c 

1+K-t Pp j-X P 

■ IlnTj’ T+^ J = - X » ^ ^ 


-b 

f 

-2 a 
< -b) , 


where it is assumed that the stiffeners are located on {-2a < x < -b, y=0) 
and (b < X < 2a, y=0), and the crack on (x=0, -d < y < -c). The equilibrium 
and the single- valuedness conditions under which (54) and (55) must be 
solved are 



-b 

► 

f3(t)dt = 0 . 
-2a 


(56a,b) 


In this problem, too, the interesting case is that of c=0, b=0 
corresponding to the crack initiation and growth from a broken cover 
plate. In the case of a single cover plate -2a < x < 2a, in the absence 
of any cracks in the elastic half plane, the maximum tensile stress 
in the cover plate would be at x=0. If the cover plate fails at 
this point, then the problem reduces to that discussed in [3], where 
it was shown that the stress state (in the half plane) at the point 
(x=0, y=0) has a strong singularity. This would greatly enhance the 
possibility of crack initiation in the elastic half plane around this 
point. For b=0 and c=0 the asymptotic behavior of the solution of (54) 
and (55) may again be examined by defining f^ and f 3 as in (36) and (38) 
and by using the function -theoretic method described in the previous 
section. In this case, the characteristic equations are found to be 
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COtiro) = 0 , 


COtna = 0 , 


(57a, b) 


~ 


<^‘5("2a) 

(2a 


1 4*1 (O) o 4*‘3(Q) 

[(l-COSirB-46+232) + 2 (b-2)cos f = 0 , (58) 

iTmJs <=“ T — - (1+cos.e) ] = 0 . (59) 

Equation (57) gives again the standard result u = 1/2, a = 1/2 and the 
coefficient determinant a of (.58) and (59) may be shown to be 

A(6) = . (60) 

sin^irg 

It is again seen that a(b) = 0 has no root (for which 0 < Re(e) < 1). 

If one now assumes that at the end point t=0 f-j and f2 are bounded 
and are of the form given by (49a) and (49c), by using the procedure of 
the previous section, it may be shown that 62(0) and hence f3(0) is 
zero. In the asymptotic relation obtained from the second integral 
equation, the coefficients of the logarithmic terms cancel out, imply- 
ing that ((>1(0) is bounded and (may be) nonzero (see, for example, [9] 
for details). 


7. Results. 

The integral equations with a simple Cauchy kernel or with a gen- 
eralized Cauchy kernel found in the previous sections are solved by 
using the numerical integration formulas described, for example, in 
[7]. The stress intensity factors defined by (27-29) for the problem 
of a uniformly loaded stiffened cracked elastic half plane are given 
in Tables 1-6. Tables 1 and 2 show the stress intensity factors for 
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the interface shear stress at the end points x=0 and x=-2a of the stiff- 
ener for the case of an internal crack (c > 0, e = 0, Figure 1), In 
Table 1, the length of the crack is fixed and its relative distance to 
the boundary s^ = (d+c)/2a is varied. In Table 2, the crack length is 

varied for a fixed distance. In the tables the dimensionless constant 
X* defined by 

^ (61) 

is a measure of the relative stiffness of the cover plate, smaller x* 
corresponding to stiffer cover plates. The tables indicate the expected 
trends, namely that generally the stress intensity factors increase 
with increasing cover plate stiffness. Tables 3 and 4 show the corre- 
sponding stress intensity factors at the ends of the internal crack 
which are defined by (27 and 28). It may be seen that the mode II stress 
intensity factors k 2 (c) and k 2 (d) are very small in comparison with the 
mode I values k^(c) and k^(d) indicating that a subcritically growing 
crack generally would remain in the direction approximately perpendicular 
to the boundary. One may also observe that the mode I values tend to 
slightly decrease with increasing cover plate stiffness or decreasing x*. 

Table 5 shows the results for the edge crack (i.e., c=0) in a half 
plane stiffened by a single cover plate. Note that in this case, k2(0) 
for the interface shear is zero. The results are given only for (d/a) < 

1 as they appear to remain relatively constant for (d/a) > 1, This may” 
be seen from Figure 5 giving the mode I stress intensity factor at the 
crack tip as a function of d/a. The asymptotic value of the stress 
intensity ratio for (d/a) -»■ « shown in the figure is that of a uniformly 
loaded unstiffened half plane containing an edge crack. The figure 
shows that as (d/a) 0 the stress intensity factor becomes unbounded. 

This is, of course, due to the fact that in this case, the governing 
stress field itself is singular (see Figure 3). It may also be seen 
that as X* decreases, the stress intensity factors tend to increase due 
to the increase in the "stress concentration" around the stiffener. 
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The results for two symmetric cover plates are shown in Table 6. 

In this case, the stress intensity factor for the contact shear kgCO) 
and that at the crack tip for mode 11 kgCd) are zero. The important 
stress intensity factor k^(d) is also given in Figure 6. Again, it may 
be seen that for (d/a) > 1, the effect of the cover plates appears to be 
negligible (this general result may also be observed from the stress 
distributions given in Figures 3 and 4). Comparison of the results 
given in Figures 5 and 6 indicates that the stress intensity factor for 
the two cover plate case is consistently greater than that for a single 
cover plate. 
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Table 1. Variation of the stress intensity factors 

l< 2 (- 2 a) and l< 2 ( 0 ) for the interface shear stress 

with the crack location, SQ=(d+c)/(2a) - the case 

of internal crack, X*=(a/h)(y«/y,)(l+K,)/(l+K«), 
£=(d-c)/(2a) = 1. ^ 2 



Table 2. Variation of k2(-2a) and k2(0) with crack 
length, £.= (d-c)/2a, SQ=(d+c)/(2a) = 1. 
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Table 3. Variation of the stress intensity factors at the crack tips s=c and s=d with 
the crack location, SQ=(c+d)/(2a), A=(d-c)/(2a)=l (internal crack) , 



I 

ro 

Ca> 

I 


Table 4. Variation of the stress intensity factors at the crack tips with the crack 
length, Ji=(d-c)/(2a),SQ=(c+d)/(2a)=l (internal crack). 



ki(c)/pQ/iia 1 

k2(cJ/Po/iia 


K2va;/Po*'*a 

A 

0.2 

1 

4 

0.2 

1 

4 

0.2 

1 

4 

0.2 

1 

4 

0.1 

0.913 

0.935 

0.967 

0.035 

0.027 

0.015 

0.919 

0.940 

0.970 

0.039 

0.030 

0.016 

0.25 

0.919 

0.943 

0.979 

0.033 

0.026 

0.014 

0.932 

0.952 

0.982 

0.042 

0.032 

0.017 

0.5 

0.965 

0.993 

1.037 

0.023 

0.020 

0.014 

0.966 

0.987 

1.018 

0.050 

0.039 

0.021 

0.75 

1.130 

1.160 

1.212 

-0.008 

0.001 

0.010 

1.021 

1.044 

1.081 

0.062 

0.050 

0.030 

0.9 

1.500 

1.512 

1.549 

-0.083 

-0.050 

-0.004 

1.080 

1.105 

1.145 

0.073 

0.061 

0.040 



Table 5. 


k2(-2a)/pQV^ 


Stress intensity factors l<2(-2a), k^(d), and k2(d) for the case 
of a single stiffener and an edge crack (c=0) (see insert in Figure 5), 
I k^(d)/pQvW2 I k2(d)/pQ/372 . 


\ A’" 

d/a\ 

0.2 

1 

10 

0.2 

1 

10 

0.2 

1 

10 

0.1 

0.441 

0.369 

0.169 

2.670 

2.434 

1.838 

-0.380 

-0.300 

-0.095 

0.25 

0.393 

0.335 

0.161 

1.959 

1.877 

1.667 

-0.184 

-0.143 

-0.039 

0.5 

0.308 

0.269 

0.138 

1.687 

1.667 

1.606 

-0.071 

-0.055 

-0.015 

1 

0.152 

0.138 

0.081 

1.591 

1.589 

1.586 

-0.010 

-0.008 

-0.002 


Table 6. Stress intensity factors k2(-2a) and k-|{d) for two symmetric 
stiffeners and an edge crack (c=0) (see insert in Figure 6). 


k2(-2a)/pQ/a 



0.2 

1 

10 

0.469 

0.386 

0.171- 

0.402 

0.341 

0.162 

0.308 

0.269 

0.139 

0.152 

0.138 

0.081 



l<l(d)/PQ/a7? 


1 











APPENDIX A 

The functions appearing in equations (1): 

y^y-[ y-y^ 2(x-x^)^Cy+yi) 

(y+y-|)^+(x-x-j)^ (y-y-|)^+(x-x^)^ C(y+y'i)^+Cx-x^)^]^ 


2(x-x^)^(y-y^) ^ 6(x-x^ )^(y>-yT )^+2y^ (y+y^ )^-6(x-x^ )2y^ (y+y^ )-2(x-x^ 

C(y-y^)^+(x-x^)^]'^ ^ C(y+yi)^+(x-x-|)2]3 

(x-X^)2-(y+y^)2 

C(y+yi)^+(x-x^)2]2 

(x X )f (x-X^)2-(y+y^)2 

[(y-y,)2+(x-x,)2]2 ’ [(y+y,)2+(x-x,)2]2 



y[(y+yi)^-(x-x^) 2 ] + 2(y+y^)[(x-x^)2-(y+y^)2+2y^(y+y^)] 

C(y+y-i)^+(x-x^)2]3 


} 


9 


x-x-j ^ x-x-j 2(x“X-j )(y-yi )2 

(y+y^)^+(x-x^)2 (y-y^)2+(x-x^)2 C(y-yi) 2 +(x-x^) 2]2 

2(x-x^ )[(y+y^ )2-2y^ (y+y^ )] 

C(y+y-,)^+(x-x^)2]2 


(x-x^ )3 (2y+y^ )+(x-x^ )[3y^ (y+y^ )^-2(y+y^ )3] 

[(y+yT)2+(x-x|)2]3 



(y-y-i)C(y-yi)^-(x-x^)2] (y+yi)[(y+yi)^-(x-xi)2] 

[(y-yi)Mx-x,)2]2 C(y+y,)2+Cx-x,)2]2 


2yi 


[(y+yi)^+(x-xp2]3 


{4(x-x^ )2[(y+y^ )2-(x-x-j )^-2y^ (y+y-, )] 


+ [(y+y-j )^+(x-x-j )^][3(x-x^ )2-(y+y^ )2+2y-| (y+y-j )]} 


K. 


31 


K 


32 


[(y-y,)2+(x-x,)2]2 

[(y+y-l)2+(x-x^)2]2 

2yi[( y+yi)2-(x-x^)2] 

4yy^ (y+y^ )[(y+y^ )2-3 (x-x^ )2] 

[(y+yi)^+(x-x^)2]2 


C(y+y,)2+(x-x,)2]^ 

X-X^ 

X-X-J 

2(x-x,)(y-y,)2 

(y-yi)2+(x-x-,)2 (y+y-|)2+(x 

-x,)2 [{y-yjj^+Cx-x,)^]^ 


2(x-x^ ) (y+y^ 4(x-x^ )yy^ [3(y+y^ )^-(x-x^ )2] 

[(y+yi)2+(x-xi)2]2 C(y+yT)2+(x-x^)2]3 
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APPENDIX B 


The kernels k.. appearing in equations (12): 

■ J 

kyiCsst) = + [ 2 (t-s)sin 20 cose - (t+s)cos3e]R'^ 

+ [(t-s)(5t2-s2)sin20cos30-2t(t+s)2cos^0 - (t-s)3cos0sin‘*e]R"^ 
+ [4t^(t+s)3cos^0 + 8t2(t+s)(4ts-t2-s^)sin20cos20 
-12t^(t-s)2(t+s)sin^0cos30]R"^}, 

2y 

l<12(s.t) = (t-s)sin30R ^ 

+ sin0cos‘*0[(t+s)(t2-s2-8ts) + 4ts(t-s)]R“'^ 

+ sin30cos^0(t-s)(17t^+s2-lOts)R~^ 

+ [8sin0cos®0t(t-s)2(t+s)2 - sin^0cos204t(t-sf(4t-3s) 

+ 4sin30cos^0t(t^-s^) (s2-2t2+7ts)]R~®} , 

O o 

kio(s,x^) = - -H[(s sin0-x^) cos^0 + s^sin20cos^0(s sin0-x.) 

13 0 TT*- 0 0 

-2sin0cos20s(s sin0-XQ)^][s2cos^0+(s sino-x^)^] 
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2y 

^ ^(<2+1 ) J^^t-s)(.s1n39-sin0cos2e)-Ct+s)s7necos20]R"2 


+ CCt+s)C2ts-s2-5t2)sin0cos‘*0 + (t-s)2(3t-s)sin30cos2e]R'^ 

+ [4t(t+s)2(t-s)C2t-s)sin0cos60 + 8ts(t-s)(3t2-s2)sin30cos'*0 

-4t(t-s)3(2t+s)sin^9cos20]R“®J , 

2y 2 

' ^^T+Kg)' " (t-s)sin20cos0R-^ + (8t(t-s)2cos0sin‘+0 

+ cos30sin20[(t+s)(t2-s2-4ts) - 2(t-s)(5t2+s2-2ts)] 

+ (t+s)(s2-t2+4ts)cOs50)R“‘* + [4t(t-s) (t+s) (4t2-3s2 

- 7ts)sin20cos50 + 4t(t-s)3(2t-s)sin‘*0cos30 

“ 8 t(t-s)‘+sin 60 cos 0 ]R"® } 

9 

k23(s»Xo) = - I [s2cos20 + (s sin0-XQ)2]"2[(s sin0-XQ)^sin0cos0 

- s2(s sin0-XQ)cos30sin0 - s(s sin0-XQ)2(sin20-cos20)cos0], 

R2 = (t+s)2cos20 + (t-s)2sin20 . 
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k 3 i(x,t) = {-T“^ + [t2cos2e - 3(x-tsine)2]T’'*}, 


l< 32 (x,t) = - j t2(x-tsine)cos2e , 


SHil 


H(x-x„) , 


T = t^cos^e + (x-t sine)2 , 


1 , x > x_ 
H(x-x ) = 1 ° 

0 , X < Xq • 
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Fig. 1 



Geometry of the stiffened elastic half plane containing 
a crack. 




0 


Q 


7T/2 


^/2V/ 


7T/4 


tt/4 


Fig. 2 Angular distribution of the cleavage stress in the half 
plane around the end point of the stiffener. 





Fig. 4 Variation of the normal stress on the surface of the 

aX 

stiffened half plane. 
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Fig. 5 Variation of the mode I crack tip stress intensity factor with 
the crack length for various stiffness parameters x* = 
[ay2(l+<i)]/[hy-i(l+K2)3 in a half plane containing an edge 
crack and a single stiffener. 







